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Abstract 



The Schwinger model is used to study the artifacts of quenching in a con- 
trolled way. The model is solved on a finite-temperature cylinder of circum- 
ference j3 = l/T with bag-inspired local boundary conditions at the two ends 
x 1 = and x 1 = L which break the 75-invariance and thus play the role of a 
small quark mass. The quenched chiral condensate is found to diverge expo- 



^^ ■ nentially as L — ► 00, and to diverge (rather than melt as for JVf > 1) if the 



high-temperature limit /3— >0 is taken at finite box-length L. We comment on 
the generalization of our results to the massive quenched theory, arguing that 
the condensate is finite as L — > 00 and proportional to l/m up to logarithms. 
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I. INTRODUCTION 

Simulations of lattice QCD often use the quenched approximation in which the quark de- 
terminant is replaced by a constant []]]. While it is clear that the quenched theory is 
unphysical 0, it is less certain what pathologies are introduced in Euclidean correlation 
functions. Considerations based on quenched chiral perturbation theory suggest that the 
chiral limit is singular ||. For example, the quark condensate is predicted to diverge in 
this limit: (ipip) oc m~ s , with 5 ~ 0.1 []. There is some numerical evidence supporting the 
predictions of quenched chiral perturbation theory, but it is not yet definitive || . 

Given this situation, it is interesting to explore the effect of quenching in a completely 
controlled environment where calculations can be done analytically. In this paper we un- 
dertake such an investigation for the Schwinger model (QED in two dimensions with Nf 
massless fermions ||). This has been frequently used as a toy model for QCD(4), since it 
shares the properties of confinement and dynamical mass generation J7J. In particular, we 
have calculated the chiral condensate in the massless quenched theory with chiral symmetry 
breaking introduced by the spatial boundary conditions. The chirality violating parameter 
is 1/L (L being the box length) instead of the quark mass m. We also introduce thermal 
boundary conditions in the Euclidean time direction. In this way we can formally work with 
the massless theory throughout (the finite spatial and temporal extent of the manifold act 
as an infrared regulator), and maintain full analytic control. In particular, the result for 
the theory with positive N{ quantized with these boundary conditions can be analytically 
continued to the quenched limit, N{ = f\. We are able to calculate the condensate in the 
quenched theory at any spatial location as a function of L and the inverse temperature (3. 

A more straightforward approach would be to quantize the massive theory with an ar- 
bitrary number of flavors in a finite box, and then take the limits N{ — > 0, L —* oo and 
m —* sequentially (in this particular order). The difficulty here is that even in the case of 
the massive QED (2) only an approximate solution exists ||. Since analytical continuation 
in iVf is, in general, arbitrarily sensitive to small variations in input data, one is on more 
secure ground using an exact solution. The drawback in our approach of breaking the chiral 
symmetry by boundary conditions is that the last two limits mentioned above (L — > oo and 
m — > 0) are effectively being taken simultaneously. 

There have been several previous investigations of the quenched Schwinger model, both 
analytical and numerical. Casher and Neuberger argue that the spectral density of the 
Dirac operator at zero eigenvalue is non-zero || - - an observation which, if we were in four 
dimensions, would indicate spontaneous chiral symmetry breaking. Introducing an ad-hoc 
infrared regulator mass /irR, Carson and Kenway [|K| and Grandou |[l| conclude, using 
bosonization methods, that the quenched condensate actually diverges exponentially as the 
regulator is removed: 



1 This divergence is not related to that due to exact zero modes in topologically non-trivial back- 
grounds, an effect which is suppressed in the infinite volume limit Q. 

2 In the presence of an infrared regulator, this is equivalent to directly calculating correlation 
functions in the quenched approximation, i.e. without the determinant. 
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(i/j^i/j) ~ /x m exp - — T . (1) 

Smilga reaches a similar conclusion based on an analysis of the eigenvalue spectrum of the 
Dirac operator [[E|. Our calculation shows that the result ([[]) is just a specific form of the 
generic combination "power-law fall-off times exponential divergence" , different versions of 
which we find as exact asymptotic results when the boundaries are sent to infinity in various 
ways. Our approach uses a finite box with specific chirality breaking boundary conditions 
rather than an ad-hoc IR-regulator and allows for a generalization to finite temperatures. 

The following section summarizes previous work on the Schwinger model with chirality 
breaking boundary conditions JT3-irS[. We then describe the analytic continuation to N{ = 0, 



and study its behavior in the interesting limits. In sec. [TV| we elucidate the origin of the IR 
divergences that we find. Section |V| is devoted to a more speculative discussion of the IR 
divergences in the massive quenched theory. We conclude with a summary of our results. 

II. CHIRAL CONDENSATE WITH ARBITRARY NUMBER OF FLAVORS 
A. Manifold Parameters and Abbreviations 

The Euclidean Schwinger model (massless QED in d = 2 dimensions) 

S[A,^^]=S B [A] + S F [A,^^] 
1 N f 

- A JF» V F, U , S F = Y,JHP1>« (2) 



1 ' 
Sb 

M n=i M 



n=l >' 



is studied on the manifold 

M=[0,/3] x [Q,L] 3 (xV 1 ) (3) 

with volume V = f3L. In Euclidean time direction, the fields A and ip are periodic and 
antiperiodic respectively with period (5. Hence x° = and x° = (3 are identified (up to 
a sign) and the manifold is a cylinder. At the two spatial ends of the cylinder (i.e. at 
x 1 = and x 1 = L) some specific chirality-breaking (XB-) boundary-conditions (which will 
be discussed below) are imposed. We use the one-flavor Schwinger mass 

" : - 4 (4) 

to simplify our notation, and introduce the dimensionless inverse temperature and box-length 

a = fi/3 A = jiL (5) 

as well as the dimensionless volume, aspect ratio, and spatial position 

v = a\ , r = f x , e = Y"- (6) 



B. Chirality Breaking Boundary Conditions 

The proposal to study both QCD and the Schwinger model with chirality-breaking 



boundary-conditions goes back to Ref. ||16|| . The XB-boundary conditions can be motivated 
by requiring the operator ilj) to be symmetric under the scalar product (%, ip) :— J % y d 2 x, 
which leads to the condition that the surface integral § y}iji ip ds vanishes, where i/i = n fl , ~f ll 
with n^ denoting the outward oriented normal on the boundary. Imposing local linear 
boundary conditions which ensure this requirement amounts to having xV"^ = on the 
boundary for each pair. A sufficient condition to guarantee this is to require all modes to 
obey ip = Bip on the boundary, where the boundary operator B (which is understood to 
act as the identity in flavor space) has to satisfy B^t/iB = —rji and B 2 = 1. In ||i3|--|l5|j the 
one-parameter family of solutions 

B = B e := i^e^i (7) 

was chosen, supplemented by suitable boundary conditions for the gauge-field. The 75 
invariance of the theory is broken for all 9, thus making the iVf-flavor theory invariant under 
SU(Nf) v instead of SU(Nf) L x SU(Nf) R . In physical terms, these boundary conditions 
prevent the t/(l)-current from leaking through the boundary, since j ■ n = ip^ijiif) = on 
dM. There are considerable differences regarding the spectrum of the Dirac operator in the 
theory with XB-boundary-conditions |TB,rj] as compared to the theory on the torus ||17|j : 



- The Dirac operator has a discrete real spectrum which is asymmetric with respect to 
zero. 

- The spectrum is empty at zero, i.e. the Dirac operator has no zero modes. 

- The instanton number q = e/(4Tc) ■ J e^F^ = e/(2ii) ■ J E G R is not quantized. 

The first property already indicates that we are away from the usual Atiyah-Patodi-Singer 
index-theorem situation. The fact that the second property is fulfilled implies that the 
generating functional for the fermions in a given gauge-field background 

Z F [A,r}\rj\ = ^ / 1] 5^^%) e -S/%)W«+E/^)W+E J%)*w ( 8 ) 

does indeed simplify to the textbook formula 

from which the one-flavor condensate in a given background [no sum over (i)] 

(4 ) (x)P ± ^ ) (x)) = -L f t Z F I (10) 



is computed. Here and in the following P± = 1/2 • (1 ± 75) denotes the projector on the two 
chiralities, where 7s = diag(l, —1) in the chiral representation of the Dirac Clifford algebra. 



C. General Result at Arbitrary Points 

The condensate in the iVf-flavor Schwinger Model on a finite-temperature cylinder with the 
XB-boundary-conditions (|7|) at the two spatial ends was found to read [fL4] , |15|1 

(lp j P±?P) e ±e-cosh(A v / 7Vf(l-25)/2)/cosh(A v / 7Vf/2) 

fi ~ 4A 

1/2 

J dc cos(27mc) 9 3 f (c, it) 
sin(7r£)cosh(7mr) -1/2 



U ' sin^fl+sinhW) f ^ 

-1/2 

exp {^ E (l ~ co S ( 27 rne)) ( ° oth ^ - ( w -» ^+ iV f (A/tt) 2 )) } (11) 

iV f n>l " 

(lp i! P±1p) e ±e-cosh(A v / 7Vf(l-25)/2)/cosh(A v / 7Vf/2) 

// 2(X 

1/2 

/ dccosh(27r(m+Oc/r)e- Jv ^ c2 / T ^ t (ic/r,z/r) 



sinh(7r(m+0/r) V 2 



/ dc e- N ^ c2 /-e^(ic/r,t/r) 
1/2 



7T /£(!-£) cosh(Ax/iVf(l-20) - cosh(A N /iVf 



f 7T / ^-<U cosiHAyiVfii-z^j -cosn^yiVfj \-| 

6XP W f V r av^Vf sinh(AViVf) " ' l J 

r 1 ^ cosh( irm / r) — cosh( irm( 1—2^)/ t) , 1 — - — — — .„ NnN i 

exp F E " • u , \ JL ~ L ~ ™ ~* V m + ^f(a/2vr)2 )} , 

-^i m>\ msvahyKm/T) v ' 

where the integration variable c represents the harmonic piece in the Hodge decomposition 
of the gauge potential, and the elliptic function 9 3 is defined as 

3 (u, u) = Y, e 2 ™V 2 = 1 + 2 E cos(2nvTM)g n2 (g = e i7rLJ ) (13) 

ngZ n>l 



for the parameters uj = it and a; = z'/r (giving real nome g e]0, 1[). The two forms (|TT, [12]) 
are identical for any finite <r and A, but they enjoy excellent convergence properties in the 
regimes r ^> 1 and r<l, respectively (r = <r/2A). 

III. QUENCHED CHIRAL CONDENSATE 

A. Quenched Chiral Condensate at Arbitrary Points 

This paper is based on the observation that - even though certain expressions within formulas 
( |TT1 , |T2|) diverge as N { — > (for fixed a, A) - the product of all factors stays finite: 

(i[>iP±i[>) ^e ±e 1 ( \ 2 , ..wcoth(Trnr) 



/' 



e 1 r A v-, / , „. \ /coth(7mr) r \-. . ,. 

7Y ^ T -r-exp{— Efl-cosSTrnOjf hr^ + ^T2? 0/ 14 

4A sm(7r£) ^ ri>i nn n 2 smh (7mr) /J 



e ±e t r rca 2 



-p{^ 2 d-« 2 } 



/i 2a sin(7r£) ^ *■ 12r 

eX P IzZZ L 2 ■ ,2/ ZZl 1 - £ COsh ( ) - (W)cosh(- — 

" 7rr m>i ?7i sinii {irm/T) v r r 

r (j2 v^ 1 / , / ffm N , /7"™(1-2£)n\ /,kn 

ex P nL^T/ rz( cosn ( )-cosh( )) , (15) 

"^ m>i TTi- sinJi(7rm/r) v r r 7 

where the low- and high-temperature forms fll~4] , p~5|) are identical for any finite a, A. 

B. Specialization to the Quenched Case at Midpoints 

To study chiral symmetry breaking with XB-boundary conditions, one has to determine 
whether the condensate in the bulk of the cylinder survives when the boundaries are sent 
to infinity We choose to consider the midpoint (£ = 1/2), and also to set the boundary 
parameter 9 to zero, for which the resulting expressions are most simple. Other values of 
£ and 9 lead to the same conclusions. Specializing eqs. ( |i~4| , |15"D to £ = 1/2 and 9 = 0, and 
rewriting them in terms of the variables v = aX and r = cr/2A, we find 

(jj t P±ij} _, Vr r v y> / coth(7r(2n + l)r) r U(16l 

/i 2 3 / 2 v ^7' eXP l27rr^o^ vr(2n + l) 3 (2r2 + l) 2 sinh 2 (7r(2n + l)r)^ l ' 

WP±i/>) i \/t (jn^_ vt y, / tanh(7rm/2r) 1 x-, 

/i 2 3 / 2 v^7' eXP l96r 2 4^^/ ™ 3 2m 2 rcosh 2 (7rm/2r)^ ' l ' 

These two expressions enjoy excellent convergence properties in the regimes r ^> 1 and 
t <C 1, respectively. 

C. Asymptotic Expansions for t>1 and r<l 

From (|T6| , |TTD the condensate is seen to behave asymptotically like 

{fP±f} . 1 fA 2 ^ 1 1 , 1 f 7C(3)A 2 1 low temp 

^^^ ± 4T eXP l^g(27mF} = ± 4A eXP l^ 2 -} ( fixed L } (18) 

^^-±-exp{^Ez-^}-±-exp{»} ( L ^°) (19) 

H 2a P U7r 2 r 2 ^ (2n+l) 3 J la l I 32tt 2 t 2 J V fixed T ; 

(^tp ± ^,) : r _ f vra 2 a 2 ^ 1 1 , e C(3)<T 2 /8T 2 T ^ ^ f ^2 ^ /L ^oo 

fl ^u '^' "" m > 1 

(if>ip ± if>) 1 rvrA 2 A 2 r 2 _ 1-, 1 / ttA 2 1 high temp 

±77 ex P 1 7ZZ + IZZ^ 2_. ZZ^ r~ ± 7T ex P 1 7ZZ ( c , r )( 21 ) 



r ( na 2 a 2 „ 1 , e^^V r tt^ 2 , ,L->oc, , N 

± 27 eXP tl92^ + 8^^ 1 ^} = ± 2a 6XP 1 192^1 (fixed 2^ (20) 

m>l 

1_ rvrA 2 AV^ ^| 1 rvrA 2 . 

// ^4A eXP l487 + ^ 2 "^ 1 ^ 3 "J 4A 6XP ^487J ' fixed L 

where only the leading term is kept. 



D. Divergence of the quenched chiral condensate in the chiral limit 

As noted in the introduction, we are particularly interested in the behavior of the condensate 
in the chiral limit. The condensate is properly defined by first taking the infinite volume 
limit and then sending the chirality breaking source to zero. In our framework these two 
limits are taken simultaneously by sending A — > oo. This means that in general we do not 
know a priori whether the source is "strong enough" to lead to symmetry breaking when 
A — > oo. For example, in QCD(4) if m — > and L — > oo such that mV — > in the limit, 
with V the Euclidean volume, then chiral symmetry breaking will not occur. But if we find 
a non-zero limiting value for the condensate, we conclude a posteriori that the symmetry 
breaking source is strong enough. Indeed, we shall see that in the (massless) quenched 
Schwinger model the condensate diverges as A — > oo for any temperature. 

1. To study the limit of infinite box- length at zero temperature we first perform the limit 
of zero temperature in (|T8|) which gives the exact result 



lim 



<^> =± i { m.}. ,22, 

//. A fixed 4A •- 8vT 2 > 



Clearly the infinite volume limit cannot be taken. 

2. To study the limit of infinite box- length at fixed (non-zero, non- infinite) temperature 
we use (p0|) from which we see that the condensate now diverges like 

(±ML eC^ {JL } A»l 

y, 4A ^24a J V fixed ; v ; 

which is even more virulent than in the zero-temperature case (J2"2"D. 

3. To study the limit of infinite box- length at infinite temperature we first focus on the 
high-temperature limit at finite box-length. From fl2"Tp we find the result 

(i 4A P l 24 a J V A fixed 7 V ; 

which diverges at infinite temperature even at finite box-length. 

4. Finally one may check whether the results are significantly altered if the limit of infinite 
box-length is taken while simultaneously lowering or raising the temperature in such 
a way that either the aspect ratio r = o~/2A or the box-volume v = aX stays constant. 
Under v — > oo at fixed r the condensate is seen from fll6|) to diverge like 

(iplP±ip) 1 r i 

~ ±const —= exp \ const v \ (25) 

Under r — > at fixed v the condensate is seen from (|TT| ) to diverge like 



Our findings are illustrated in figure [l|. 



±const y/r exp I const —> . (26) 



FIGURES 




FIG. 1. The quenched condensate at the midpoint as a function of spatial box size and inverse 
temperature, using eqs. ( p5| , p6| ). The axes have been rescaled so that the entire range of values 
is shown: x = 2/ir ■ arctan(cr) and y = 2/ir ■ arctan(A), where a = fi[3 and A = fiL. Unlike the 
other divergencies, the one near y = is unphysical as it is an artifact of the boundary conditions. 
The border of the surface at x = 1 represents the (massless) quenched Schwinger model at zero 
temperature with our boundary conditions. 



IV. ORIGIN OF INFRARED DIVERGENCES 



It is interesting to elucidate the origin of the IR divergences in the quenched theory. 
Mathematically, things can be traced back to the fact that in two dimensions the gauge- 
field can be decomposed into a coexact and an exact piece (A^ = —e^ v d u (p + d^x plus, with 
certain boundary conditions, a topological and a harmonic contribution) and the functional 
determinant (if present) can be computed by integrating the anomaly. As a consequence, the 
full effective action in the iVf-flavor Schwinger model is found to read (up to contributions 
from the topological and harmonic parts) 

S nM = ^2 /^ ~ Nf^A)<p dx . (27) 

Going back to (|| || |10D and using decomposition properties of the massless Dirac operator in 
a given gauge- field (which are specific to two dimensions), the condensate with our boundary 
conditions and for arbitrary iVf > is found to factorize 



<*t(,)P ± , ( ,)) = IdcD ^. s : ix ' x) .t::2"' ,c 



fdcD<t> e~ l0 ' N f [ 
J dc S e (x, x)±± e~ N f r ^ J D<\> g^M-r^l 



Jdc e-^c) j D(j) e -v e>N} [ 



(28) 



Here Sq(x,x)±± denotes the diagonal entries of the Green function of i<j) (obeying XB- 
boundary-conditions) and Fg^Acf)] is just Sjv f [0], up to a term which vanishes for 9 = (for 
further details the reader is referred to JT^IT^]). The first integral in Q2"S| ) (which is over a 
c-number) is trivial in the quenched theory. The second integral (which is an integration 
over all 0-fields which are periodic over /3 and satisfy Dirichlet boundary conditions at x l =0 
and x l = L) can be done exactly for all Nf. 

5 r ,. =expJ2e 2 (d j— \\x) \ . 29 

f D(f> e- SN f [4>] X "A 2 -A f ^A" 7/ V ; 

We stress that the factorization property (|28|) reflects itself in the form of our analytical 
results — both in the dynamical case [eqns. (JTTJ, |TJ)] and the quenched case [eqns. (ITJL |T5|)]: 
All expressions consist of a prefactor [the first two lines in eqns. flTT| , p~2| ) or the factor 
l/(4Asin(7r£)) in eqns. (0, $Bj) ] and a product of exponentials. In both cases the prefactor 
corresponds to the first factor in (ESI), whereas the exponentials represent the second factor 



in (|28|) , as expressed in (|29j). 

The important point here is that in the quenched case the IR divergences are due entirely 
to the second factor in (pS|), i.e. due to the divergence of (^9|). Indeed, as our explicit 
expressions show, the first factor in ( PS|) vanishes as a power-law in the limit of infinite 
box-length in both quenched and full theories. This is as expected, since the leading effect 
of the chirally asymmetric boundary conditions on the free fermion propagator, S(x,x), is 
proportional to the propagator from x to the boundary and back, and this propagator falls 
off since the Dirac operator has a mass-dimension lower than the dimension of the manifold. 
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By contrast, the operator A 2 — iVf/i 2 A, which appears in the expression (p9[), has a 
mass-dimension higher than the dimension of the manifold. Its propagator 

G(x-y;N f ) = (x|| A2 _^ 2A ||y) , (30) 

thus grows with distance, and one has the possibility that the boundary conditions become 
more influential as the spatial length is increased. To understand this in more detail, it is 
useful to have explicit expressions for the requisite Green functions. In the quenched theory 
one finds a quadratic increase 

Q{x -y;N f = 0) = ^-(x-y) 2 ]og(ji\x-y\) - ^(x-y) 2 , (31) 

while for Nf > 1 the result increases only logarithmically 

G(x - y; N f > 1) = -^ ( - ^- \og(Nl /2 ^\x-y\) - ^-K (N} /2 fi\x-y\) ) . (32) 

The less rapid increase in the unquenched case turns out, as shown by the explicit results, 
to be insufficient to lead to an IR divergence in the condensate. 

One can push the argument a little further by imagining using images to enforce the 
boundary conditions. One cannot actually repeat our calculation in this way, because the 
sum over the two-fold infinite set of images does not converge, since contributions grow as 
the images become more distant - - both in quenched and full theories. But one can use 
images to obtain the pattern of IR divergences observed in our results by considering the 
simpler geometry of a semi-infinite spatial extent at zero temperature: at x 1 = we apply 
XB-boundary conditions whereas for x l — > oo we require L 2 -integrability. The analogous 
quantity to that which we have calculated is the condensate at a distance L/2 from the left 
boundary, and the important issue is how it behaves as L — ■> oo. In this new geometry we 
need only a single image, and it is the Green function between the original position and its 
image which dominates for large L. The asymptotic forms are thus 

fo 2/ „ 1 || Q Jexp{const-A 2 log(A)} (N f = 0) 

exr><2e(x\\ 3 — \\x) > = < r -, 33) 

F l "A 2 -JV f fA" /J \exp {const -log(A)} (N f > 1) 

where we use A = L\e\j ^pn instead of L to aid the comparison with our previous results. 
We expect that our results should have a similar form, but with different constants, due to 
the other mirror copies present in our geometry. 

This expectation turns out to be correct, up to logarithmic accuracy. For the quenched 
case, Nf = 0, our exact result at zero temperature, eqn. $3[ ), takes the form exp(0.107A 2 ). 
It thus matches eqn. (0), except that the logarithm is absent. For Nf — 1 the form ( |33|) 
works with "const" equal to 1, which is exactly what is needed to cancel the asymptotic 1/L- 
behavior of the first factor in (EHf) and to reproduce the well-known value for the condensate 



of the one-flavor Schwinger model in the limit L— >oo [ 14]. For Nf — 2 the "const" is found 
to be 1/2, which is not sufficient to cancel the asymptotic 1/L-behavior of the first factor in 
(HI). Consequently the two-flavor condensate vanishes as 1/a/L when the model is quantized 
with XB-boundary conditions [|15[ . 



This completes our attempt to make sense out of the functional form of our analytical 
results. While is interesting to see that the IR-sickness of the quenched Schwinger model 
can be traced back to the long-distance behavior of the Green function (x\ A~ 2 \y), it is clear 
that the argument is specific to two dimensions. 

V. SPECULATIONS ON THE MASSIVE THEORY 

Having an analytical result for the (massless) quenched Schwinger model it is interesting 
to consider the IR behavior of the corresponding massive model, quenched QED(2). In 
particular, we would like to know what happens if one disentangles the limits L — > oo, m —>■ 
and takes them in the standard order: first L — > oo and then m — > 0. In QCD(4) these two 
limits do not commute: the standard order leads to a non-zero condensate, whereas (ipip) 
vanishes when m — > is taken at any finite L. In quenched QED(2) knowledge concerning 
these two situations is asymmetric. Taking the massless limit first, analytic results are 
obtained and, as we have seen, one cannot take the L — ► oo limit when m — 0. The result 
of taking the limits in the standard order is, however, controversial. There is no exact 
analytical expression for the condensate for m ^ 0, and so one must resort to approximate 



methods. Reference [11 1 argues that sending L — > oo at fixed m yields a finite result for the 



condensate, and that this result remains finite as m — > 0. This is what is found explicitly in 



the approximate calculation of Ref. | I8| . On the other hand Ref. \T§\ proposes a power-law 
divergence as L — > oo. In the following we analyze the situation using two approaches. First, 
we expand in powers of the quark mass, and find that each coefficient diverges exponentially 



as the volume is sent to infinity. Second, we generalize the approach of Smilga ||12| , and 



argue that there is in fact no IR divergence at finite m, and that the IR divergences in our 
first approach are artifacts of expanding about m = 0. The IR divergence does, however, 



reappear as a power law divergence when m — > 0, in disagreement with Refs. f[T],|T8 . 

The derivatives of the condensate with respect to m at m = are higher correlation 
functions in the massless theory, and can be computed analytically. Note that, since our 
boundary conditions break chiral symmetry, there is no reason to expect non-analyticity 
about m = for finite volume. 

We start from the form of the condensate in the massive quenched theory 



//tn ,w ,_fdcD<p (S e , m (x,x) ++ -S ejrn (x,x) ++ + S efi (x,x) ++ ) e ^ J ^ A ^ 

( * +i ' ){) " jdcD<p e-i/« 3 * ' (3) 

where Sg tTn is the inverse of p+m with XB boundary conditions, and S$ >rn its counterpart at 
zero gauge-coupling. The subtractions are chosen such that (|34]) vanishes in the limit e — >• 



at fixed m, while still reducing to ( p8[) when m — > at fixed e. In addition, they make the 
expression UV-finite. To get a sense of the general term in the expansion about m = 0, we 
consider the first derivative, 



d(^P+ip)(x) 



dm 



rre=0 



2*12 J d V JdcSo(x,y)+±Se(y,x)±+ x 

// D(f> e - 2 <f>( x )T24>(y)-^J<i>& 2 ct> 
{ / D0 e-^I^ 2 * 
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= 2^2 j dy j dc S e (x,y) +± S e (y,x) ±+ x (35) 

(exp{2e 2 ((x\^\x)±(x\^\y}±(y\^\x} + (y\^\y))}-l\ . 

Here y runs over the entire manifold, the c-integration is the average over a c-number valued 
harmonic part of the gauge field (necessary since the cylinder is not simply connected [[I |]), 
and the sum consists of two terms (i.e. the ± are either all + or all — ). The explicit forms 
of Sq (the free massless propagator subject to the XB-boundary-conditions) and (x\A~ 2 \y) 
are given in \Ti\. What is important here is that, for geometrically fixed x and y, (x\A~ 2 \y) 



diverges as If . Thus, for at least one of the choices of sign, the second factor in fl35|) diverges 
as e L . This dominates the power-law fall-off of the first factor. Thus there is at least a 
set of y- values with nonzero measure for which the combined argument of the ^-integration 
in (|3^) is IR-divergent []. Similar arguments can be made for the higher derivatives of the 
condensate. 

We conclude that each coefficient in the expansion of the condensate about m = is 
highly IR divergent. Of course, this does not necessarily imply that the quenched condensate 
diverges when V — > oo at fixed m, because the limits may not commute. For example if the 
quenched condensate had the following schematic dependence on m and V 

(^P±*P) - - X 7 77, , (36) 

m + exp(— V) 

then it would be finite as V — > oo at fixed m, while derivatives at m = are all IR divergent. 
In the following we argue that a form like eq. (^) actually holds. 

We use the line of reasoning initiated by Smilga jl2[1 , generalized to the quenched case. 
We briefly recapitulate his argument. The key result (valid for iVf>l and for iVf = 0) is 

(E ^> = 2^3 / d " x d 'y J^yf exp{2e 2 (S(0; N f ) - G(x - y; N f ))} . (37) 

This relates the sum of the inverse-squares of the eigenvalues of the Dirac operator to the 
Green function appearing in (|3~I , p2|) . As V —> oo the leading behavior of the r.h.s. of (571) 



is determined by the asymptotic form of Q. For unquenched theories one finds, using (p2|), 

(E i> « ^ + W (iV f >l) , (38) 

with missing dimensions provided by powers of e. This implies that the characteristic size 
of the lowest eigenvalues is 



3 On this point the unquenched theory differs significantly: Each Green function diverges only 
logarithmically with L, and the overall divergence of the second factor is compensated (iVf = l) or 
overwhelmed (JVf > 2) by the fall-off of the c-integral. 

4 In the remainder of this section we use A to denote eigenvalues and not the dimensionless box 
length. 
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A c oc y-W+D/O^) (jv f >i) . (39) 

One now assumes that, in the limit of infinite V, one can define a density of eigenvalues per 
unit volume, p(\). If so, then one has approximately that 

p(A c ) « -i- oc A W-D/W+D (iV f > 1) . (40) 

A C V 

Since this holds for all large V, one can read off the small A dependence of p(X). One can 
check this result by noting that it reproduces the initial result (j38|): 



(£ ^) « / AmaX V>(A)i - yA~W +1 ) oc yW+UWr . (41) 

Here we approximate finite volume by the lower cut-off on the A integral. 

This line of reasoning can be extended to include the quenched case, N{ — 0. Smilga 
notes that if one uses the quenched result for Q(x), eq. (|31~D, one finds that, asymptotically, 
the sum rule is 

(£ h oc ex V {V} (iV f = 0), (42) 

and so the characteristic size of the lowest eigenvalues is exponentially small 

A c ocexp{-y} (iV f = 0). (43) 

We stress that in both these results we control neither the constant multiplying V in the 
exponent (aside from the fact that it is proportional to e 2 ), nor subleading power dependence 
multiplying the exponentials. 

At this point we extend the results of Smilga by determining the form of p(A) in the 
quenched theory [J The same argument used above leads to 

^> K a^ ~ 1^x7) ■ < 44) 

This cannot be completely correct because it predicts that the number of eigenvalues in an 
interval, 

n(A 2 )-n(Ai)= f" 2 dX p(X) = log{log\ 2 /\ogX 1 ] , (45) 

J Ai 

diverges as Ai —>■ 0. We propose instead the convergent form (with missing dimensions 
supplied by factors of e) 

n{X) oc -: => p(X) oc — — — . (46) 

log A A (log A) ^ 



5 We acknowledge a useful discussion with A. Smilga on this point. 
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This satisfies n(X c ) oc 1/V, which is an alternative criterion for extracting the infinite volume 
eigenvalue density from the finite volume results []. In fact, the precise form of p(X) is not 
important: all that we need is that it falls off as 1/A, up to logarithmic corrections, and that 
these corrections make it integrable as A — > 0. 

We now insert this result into the expression for the condensate, approximating the effect 
of finite volume by a lower cut-off: 

/•oo 2th 

& p *'»'"L dX JFT^« x >- (47) 

There are two cases to consider. If m < X c ~ e~ v , then the IR cut-off is provided by A c 

/■OO Orr) 777 

(^P ± ^)^j x ^r°c- (m<A c ). (48) 

If m > X c then the mass provides the IR cut-off, and one has 

rm 9 poo 9to 

(tfP ± il>) „ / dX -p(X) + / — 

./exp(-V) m Jm X 6 

2[n(m) - n{X c )\ 1 
m m 

~-[l + 0(l/V)] (m>A c ), (49) 

m 

where we have dropped subleading logarithms, which are not controlled. In other words, 
the condensate starts out at zero when m = 0, as required in finite volume with no chiral 
symmetry breaking, rises linearly until m ~ A c ~ e~ v , and then falls as 1/m up to loga- 



rithms. The point we wish to stress is that eqs. (f48| , |49| ) have no IR-divergences. Sending 
the volume to infinity at finite m, one ends up with a finite value. If one then sends m — > 0, 
however, the IR-divergence returns as a 1/m divergence. 

The preceding argument applies for boundary conditions which do not break chiral sym- 
metry, and thus is not directly applicable to the massive theory with chirality breaking 
boundary conditions. Nevertheless, since a non-zero mass removes the IR divergences with 
chirally symmetric BC, it is plausible that it will also do so with chirality breaking BC If 
so, the IR divergences we found in the mass-derivatives of the condensate must be special 
to m — 0. We conclude this section by presenting a simple model which leads to this result. 

We propose that, when considering the condensate far from the chirality breaking bound- 
aries, their effects can approximately be represented by using Smilga's analysis but with a 
small, volume dependent, shift in the mass. More precisely, we suggest that one can use 
eqs. (Ha M) with 



m — > rrixB = m + X c ~m + e v . (50) 

The shift by A c is chosen so as to reproduce the form of our exact results. In particular, it 



puts us in the region where (|49| ) holds, so that 



6 In the unquenched example above this criterion gives the same result for p(X) as obtained in 
eq. ©. 
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WP±1>)XB * = -— T-T7T • (51) 

m X B m + exp(— V) 

This simple form agrees with all our results and expectations: when m = we recover 
the exponential IR divergence of our exact result; when V —* oo at fixed m there is no 
IR divergence; and each successive derivative, evaluated at m = 0, is more IR divergent. 
This shows explicitly how the perturbative argument presented above can be correct but 
completely misleading. 



VI. SUMMARY 

We have presented a study of the chiral condensate in the (massless) quenched Schwinger 
model at finite temperature and with bag-inspired spatial boundary conditions which play 
the role of a small fermion mass. Our main analytical result is given in eqs. flTJL ITH). 



Our first finding is that the (massless) quenched Schwinger model is ill-defined due to an 
infrared embarrassment: At finite box- length with the XB- boundary conditions applied at 
the two spatial ends, it is a well-defined theory, but the condensate shows a singular behavior 
when the boundaries are sent to infinity. The precise form of this divergence depends on 
the details how the volume of the manifold is sent to infinity, but the generic structure is 
"power-law fall-off times exponential divergence" . 

Our second observation is that the condensate in the (massless) quenched Schwinger 
model does not "melt" at high temperatures, instead diverging in the limit of infinite tem- 
perature. This has to be contrasted to the case Nf > 1 where the Schwinger model is known 
to show the regular behavior, i.e. a condensate which decreases when the temperature gets 
large (see e.g. [ |T7| , |H)| and references therein). 

Both predictions can be checked using numerical simulations on a sufficiently large lattice, 
and an interesting step in this direction has recently been taken in Ref. [^C] . 

From a conceptual point of view it should be emphasized that our results for the massless 
case represent exact analytic findings. They stem from a textbook-style analytical evalu- 
ation of the path integral with fields subject to the constraints imposed by the boundary 
conditions. The latter break chiral symmetry explicitly and prevent exact zero modes of 
the massless Dirac operator. This shows that the IR divergence of the condensate in the 
quenched massless theory is not tied to a zero in the spectrum of the Dirac operator on a 
certain class of gauge field configurations. 

On the other hand, our results for the massive theory are more speculative. By extending 
the approach of Smilga |T2| we argue that the quenched condensate is IR finite, and reconcile 



this with our finding that the derivatives of the condensate with respect to the quark mass 
are IR divergent at m = 0. 
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